Temperature-induced pair correlations in clusters and nuclei 
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The pair correlations in mesoscopic systems such as nm-size superconducting clusters and nuclei 
are studied at finite temperature for the canonical ensemble of fermions in model spaces with a fixed 
particle number: i) a degenerate spherical shell (strong coupling limit), ii) an equidistantly spaced 
deformed shell (weak coupling limit). It is shown that after the destruction of the pair correlations at 
T — by a strong magnetic field or rapid rotation, heating can bring them back. This phenomenon 
is a consequence of the fixed number of fermions in the canonical ensemble. 
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The pair correlations in a macroscopic superconduc- 
tor are destroyed by increasing the temperature or the 
external magnetic field. The critical field which marks 
the boundary between the superconducting and normal 
phases, is a decreasing function of the temperature T. 
The BCS theory, which is the mean field approximation 
based on the grand canonical ensemble, describes very 
accurately this regime. Applying the grand canonical 
mean field approach to rotating nuclei |l| , where the an- 
gular velocity plays the role of the magnetic field, gives 
an analogous result: The angular velocity where the pair 
correlations disappear decreases with the temperature. 
Nuclei and atomic nano-size clusters are composed of a 
small and fixed number of particles, the single-particle 
spectrum is discrete and the level spacing is comparable 
with the pair gap. Due to these facts, the fluctuations 
of the order parameter become important, which smear 
out the transition from superconducting to normal phase 
|H HI, 0, HI and lead to pronounced differences between 
a system with even and odd particle number [(|. 

In order to properly take into account these fluctua- 
tions one has to use the canonical ensemble. The most di- 
rect way is to calculate the partition function from the ex- 
act eigenvalues of the Hamiltonian, which is possible for 
some models. In this Letter we study the simple model 
of fermions occupying an isolated shell of single parti- 
cle states and interacting by a pairing force. We shall 
demonstrate that at zero temperature the magnetic field 
or the angular velocity attenuate the paring in a step- wise 
manner until it disappears completely above a critical 
value. For a mesoscopic system in the strong magnetic 
field heating may bring back the pair correlations. This 
surprising effect is a consequence of the fixed number of 
fermions in such a small system. The reduction of the 
fluctuations in particle number leads to a strong increase 
of the fluctuations of the order parameter, which con- 
stitute the pair correlations above the critical field. A 
re-entrance of pair correlations has first been discussed 
by Balian, Flocard and Veneroni 6j , who studied ensem- 
bles with either only even numbers of particles or only 
odd numbers of particles. We shall demonstrate that the 



more stringent restriction to a fixed number of particles, 
which is realized in small systems, leads to qualitatively 
different results. 

First we consider fermions in a spherical potential. 
Then the electronic states in the cluster (the spin-orbit 
coupling can be neglected) are characterized by I, the or- 
bital momentum of the shell, its z-projection A and spin 
projection er, which we label by k = {A, a}. The condi- 
tion, k > 0, used below in Eq. means that A + a > 0. 
Due to the strong spin-orbit splitting in nuclei the nuclc- 
onic states in a shell are specified by the total angular 
momentum j, (j — I + s), and its projection m = k. 
We assume that there are only pair correlations within 
the last partially filled shell. The corresponding "single- 
shell" Hamiltonian 

H = ff pair - uM, M = L Z + gS z , (1) 
H pail = -GA+A, A+ = Y. a t a t> 

k>0 

consists of the pairing interaction H pa \ r , which acts be- 
tween the valence fermions in the last shell with the 
strength G, and the "cranking term" . The z-components 
of the total orbital angular momentum and spin, which 
are the sums of all valence fermion contributions, are de- 
noted by L z and S z . The operator A + creates a fermion 
pair in the time-reversal states (k, k). In the case of clus- 
ters, the cranking term represents the interaction of the 
electrons with the magnetic-field, where we introduce the 
Larmour frequency to = ^bB and the Bohr magneton 
fiB- The gyromagnetic ratio g = 2. In the case of nu- 
clei, the cranking term takes into account the rotational 
perturbations at the angular frequency u> and g = 1. 
The eigenvalues of H pair (uj = 0) are @ 

E„ = -j(N sh -v)(n + 2-N sh -v). (2) 

were A s h is the number of particles in the shell, which 
consists of degenerated single-particle states and fi = 
41 + 2 for clusters and 2j + 1 for nuclei. The seniority v, 
which is the number of unpaired particles, is constrained 
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by < v < A'sh. We assume that v < £1/2, otherwise 
one can pass to the hole representation. Each i^-state at 
v > 1 is degenerated. The eigenvalues of the Hamiltonian 
(HJ are 



E Vt i(ui) = E v — uM Vt , 



(3) 



The index i involves additional quantum numbers of the 
instate including the spin and possible orbital momenta 
and their z-projections for clusters or the total angular 
momentum of v nucleons and its projection for nuclei. In 
the case of clusters, /isM,,^ is the total magnetic moment 
of v particles while in nuclei M v i is the angular momen- 
tum z projection of v nucleons. For simplicity, we refer 
to M as the angular momentum and omit the index i 
for the states with minimal energy at each which have 
maximal M v . 

The canonical partition function Z and the mean value 
M are given by 



Z(7» = ^exp(- 



M(T,lu) = 



rji ) 1 

TdZ 



Z duj' 



(4) 
(5) 



The evaluation uses the fact that for each v the sum over 
i can be reduced to a sum over single-particle projections. 
For nuclei, the numerical diagonalization procedure de- 
scribed in Ref. |8( is used, which permits us also to treat 
a non-degenerate shell (see below). We consider only 
the part of M which is generated by the particles near 
the Fermi surface that participate in the pair correlation. 
The contributions for the other particles can be found in 
ref. and will be discussed together with the details of 
the evaluation of the sums for clusters in an forthcoming 
extended paper [Tof . 

We introduce the "canonical" pair gap A can as a mea- 
sure of the correlation energy 

AL„(7>) = J^5>,i|A+,4|M>exp(- E "' [G) 



Z{G) 



T 



E^< *> °l A+A ^ f ' °> exp(- ^ ( S~ 0) ),(6) 



Z(G = 0) 



T 



where \v,i,0) denote the uncorrclatcd fcrmion config- 
urations in the shell. The second term subtracts the 
expectation value of the pairing interaction in an en- 
semble of uncorrelated fermions. A detailed discussion 
of theproper definition of A can was given in the re- 
view |3j. Applying the mean-field approximation and 
the grand canonical ensemble to our model, the thus in- 
troduced A can becomes the familiar BCS gap parame- 
ter A m f. However, A can must be clearly distinguished 
from A m f because it incorporates the correlations caused 
by the fluctuations of the order parameter. We take as 
energy scale E(0) = GT2/4, the quasiparticle energy at 
T = u) = 0. A value of E(0) = 0.3 - 0.4 meV was found 
for Al-clusters with radii R = 5 — 10 nm [TTj . which sets 
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FIG. 1: Canonical gap A can (T, 10) (full lines) and the mean- 
field gap A m f(T, uj) (dotted lines) v.s. the frequency to for a 
spherical shell. 
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FIG. 2: Canonical (T, a;) for even (full lines) and 

odd (the dashed line) particle number, and the mean-field 
gap A m f(T, uj) (dotted line -BCS) v.s. the temperature T for 
a spherical shell. 



the energy scale. The nuclear mass measurements give 
E(0) = 0.8 — 1.5 MeV. Let us first consider the destruc- 
tion of the pair correlations at T = 0. Since M v increases 
with v, the energy E u (u)) becomes smaller than E v ^2{<+>) 
at the frequency 



uj v (M v - Af„_ 2 ) = E v - E^ 2 - 



(7) 



When the state of the lowest energy changes from v — 2 
to two more electron states are blocked and the pair 
correlations are reduced. The last step at the critical fre- 
quency Wcrit leads to the maximal seniority ^ max , where 
all particle states are blocked and the pairing is com- 
pletely destroyed. Fig. 2] illustrates the step-wise de- 
struction of pairing by blocking for a half-filled / = 5 
shell, which is the Fermi level in a spherical Al-cluster 
with about 10 3 atoms. 

Fig. n also shows the gap A m f(w,T) obtained by ap- 
plying the mean-field approximation (cf. 0, 0) to the 
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single-shell model. For T = 0, the pair correlations are 
more rapidly destroyed than for the exact solution. The 
quantal fluctuations of the order parameter generate ad- 
ditional pair correlations. For finite T, the mean-field 
gap behaves as known from macroscopic superconduc- 
tors: A m f(cj,T) is reduced at T — 0.4E(0) and reaches 
at a lower value of uj. This is the familiar shift of the crit- 
ical frequency w cr j t toward smaller values with increasing 
T. However, the canonical gap A can (aj, T) behaves differ- 
ently: The abrupt drop around w cr i t is smoothed out by 
the fluctuations of the order parameter. Moreover sub- 
stantial pair correlations appear in the region above uj cr \t, 
which increase with T. The comparison with A m f(o;,T) 
shows that the fluctuations contribute more at finite T 
than at T = 0. For T = E(0), the pair correlations fall off 
very gradually with u>. The nuclear case is quite similar. 

Fig. [3 shows how these temperature-induced pair cor- 
relations appear with increasing T . For uj = 0, there is a 
pronounced drop of A can around T a - lt {uj = 0) = E(0)/2, 
where the mean-field gap A m f goes to zero. Above this 
temperature there is a long tail of pair correlations caused 
by the fluctuations. For w = 0.05 E(0), the drop is 
shifted to smaller T by about the same amount as the 
T cr it of the mean-field solution (not shown). This trend 
continues for ui = 0.1 E(0). For larger uj > oj cr it, the 
pair- correlations built up with increasing T. 

The temperature-induced pairing can be understood in 
the following way. At T — 0, all electrons are unpaired 
when the state of the maximum seniority becomes the 
ground state for uj > w C rit- At T > 0, excited states 
with lower seniorities enter the canonical ensemble, which 
reintroduce the pair correlations. 

The degenerate spherical shell corresponds to the 
strong coupling limit (A/d > 1, d distance between the 
levels) of ref. [fj. For uj = 0, one may compare A can 
of the present work with Ap obtained for the ensemble 
with good particle number parity in ref. which we 
will refer to as BVF in what follows. For T < 0.5 T cr i t , 
both are similar. However for T > T cr jt, Ap = whereas 
A can remains finite up much higher temperatures. The 
case uj ^ is not shown in ref. [6| for the strong coupling 
limit. 

In order to investigate the consequence of deviations 
from the spherical symmetry, we used the deformed shell- 
model described in |8|]. The Hamiltonian is given by 



ujM. 



(8) 



The solutions are found by numerical diagonalization in 
the configuration space of a j-shell. 

Figs. and 0] illustrate the case of a cluster without 
spherical symmetry. The half- filled shell consists of 
12 equidistant levels, each of which contains two states 
with spin up and down. Our choice of the level distance 
d = 0.85 E(0) corresponds to weak coupling A « d and 
is realistic for the nano-clusters. We assume an irregular 
cluster shape. As a consequence, the orbital momentum 
is quenched and only the spins contribute to the magnetic 
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FIG. 3: Canonical gap A can (T, ui) for an equidistantly spaced 
I = 5 shell in a cluster. The mean- field gap A m f(T, ui = 0) is 
shown by the dash-dotted line (BCS). 
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FIG. 4: Angular-momentum M(T, uj) for an equidistantly 
spaced 1 = 5 shell in a cluster. 



moment ^bM. The behavior of A can and M(T, uj) is sim- 
ilar to the spherical case. The step- wise destruction of 
the pair correlations corresponds to subsequent spin flips 
of electrons, which increase M by 2. The temperature- 
induced pairing appears in the region above uj ~ 2 d. The 
spin flips have been observed in tunneling experiments 
on Al nano-clusters 0- The mean- field gap A m f(T = 0) 
breaks down at the first flip, as discussed by Braun et 
al. [13 in their analysis of the tunneling spectra based 
on the mean-field approximation. However, A can (T = 0) 
is more stable and disappears only after all levels are 
blocked by subsequent spin flips. Fig. [S] demonstrates 
that the temperature dependence of A can is qualitatively 
similar to the strong coupling case in Fig. [21 A new 
feature is that sometimes A°^ > A!!™ n . For the curves 
uj = 2d in Fig. this happens because the frequency is 
above the first crossing in the even system (cf. Fig. I3J) 
but still below the first crossing in the odd system, which 
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FIG. 5: Canonical gap A can (T, ui) for an equidistantly spaced 
/ = 5 shell in a cluster with even (full line) and odd (dotted 
line) particle number. The mean- field gap A m f(T, u> = 0) is 
shown by the dash-dotted line (BCS). 

means two states are blocked in the even but only one 
in the odd system. For u> = Ad the analogous happens 
at the third crossing. In the region of high temperatures 
(T > 2T cr it) the canonical pairing gaps for odd and even 
systems practically coincide because all members of Z 
contribute with a similar weight. At small T they differ 
because only the states with small seniorities are impor- 
tant. 

Our choice of parameters lies between to the weak cou- 
pling cases A/d = 1.14 and 1.19 of the grand canonical 
ensembles with even or odd particle number studied by 
BVF (A/d — wfA in 6]). As in the case of strong 
coupling, A can decreases very gradually for T > T cr - lt , 
whereas Ap of BVF drops sharply to zero at T crit . How- 
ever, also in the region T < T cr it the results of the two ap- 
proaches are different. The differences between the even 
and odd systems are much less pronounced in the canon- 
ical ensemble. For uj = 0, A/d = 1.14, and odd particle 
number, BVF find Ap ~ for T < 0.2 T crit , A F ^ for 
0.2 T crit < T < T cr i t , and A F = for T > T clit , which 
they called the re-entrance phenomenon of pairing. As 
seen in Fig. the strong fluctuations in the canonical 
ensemble keep A can finite in the whole temperature inter- 
val, i. e. there is no re-entrance of pairing for uj = and 



odd particle number. On the other hand, the canonical 
ensemble gives temperature-induced pair correlations for 
uj > 2 d both for even and odd particle number, which has 
not been found for grand-canonical ensembles with fixed 
the particle number parity. It should be mentioned that 
BVF considered a systems with N = 100, 101 whereas 
we studied systems with N ~ 10, for which the conser- 
vation of particle number is more important. 

We have also studied the case of a half- filled j = 11/2 
shell in a deformed axial nucleus by assuming that 
in Eq. JSJ) is proportional to k 2 . We chose the distance 
e 7/2 — 65/2 = 0.28 E(Q), which is realistic. If the axis of 
rotation is parallel to the symmetry axis, the behavior is 
similar to the deformed clusters, except the steps in M 
are different from 2. If the axis of rotation is perpendicu- 
lar to the symmetry axis, the projection M of the angular 
momentum is no longer conserved. Then M(T = 0,lu) is 
no longer a step function and A can (w, T — 0) decreases in 
a gradual manner. The increase of A can (o;,T) with T is 
found to be weak. Hence, temperature-induced pairing is 
expected in nuclei that build up large angular-momentum 
by aligning the individual angular-momenta of the nu- 
cleons near the Fermi surface. These are either spherical 
nuclei (see Figs. lll2|l or the high-K isomers (see e.g. [T^). 

In summary, at very low temperature an increasing 
external magnetic field causes the magnetic moment of 
small superconducting clusters (R < 5 nm) to grow in 
a step- wise manner. Each step reduces the pair corre- 
lations until they are destroyed. However, with increas- 
ing temperature the steps are washed out and substan- 
tial pair correlations re- appear for high field strengths, 
where they are quenched at T = 0. Nuclei that built up 
angular momentum along a symmetry axis behave in an 
analogous manner. The pair correlations are destroyed 
in a step- wise manner by subsequent alignment of the an- 
gular momenta of individual nucleonic orbitals with the 
symmetry axis. These steps are washed out with increas- 
ing temperature and pair correlations appear at values 
of the rotational frequency, where they are quenched at 
T = 0. This phenomenon of temperature-induced pair- 
ing reflects the strong fluctuations of the order parameter 
in very small systems with a fixed particle number. 
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